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Simple statistical agglomeration models can provide a universal link be- 
tween the local structure and the glass transition temperature in network 
glasses. We first stress the physical features of the models and the relevancy 
of the hypothesis made and then show how to define the glass transition tem- 
perature. The models are applied to various types of binary, ternary and 
multicomponent chalcogenide glass networks and the predictions compared 
to experimental data. 
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I. INTRODUCTION 

Although much attention has been devoted to the understanding of the glass transition 
problem0, a general relationship between the temperature of this transition (when measured 
under standard conditions, at e.g. constant heating rate), and some easily reliable quan- 
tities is still lacking!. In this paper, we show that there are some aspects of structure or 
connectivity that apparently play an important part in determining the absolute magnitude 
of Tg. The construction and the prediction of this temperature from solvable agglomeration 
models is parameter-free and can be easily extended from binary to ternary, etc. glassy 
systems. 



*Proceedings of the Cargese CNRS school Physics of Glasses: structure and dynamics, held in 
Cargese, May 10 — 22*^ 1999, edited by American Institute of Physics 
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II. AGGLOMERATION MODEL 



Let us imagine a liquid that is slowly cooled and atomic motions are progressively ar- 
rested. In network glasses, the Arrhenius-like increase of viscosity upon cooling to the glass 
transition is intimaly related to a decrease of dangling bonds as the starting network is 
polymerized. One typical physical process taking place in the supercooled liquid should 
therefore be a kind of sticking process in which clusters (or macromolecules) agglomerate 
together. Also, one should remark that the most important determinant of chemical and 
physical properties of a glass is the concentration of different types of atoms involved. Thus 
the simplest level of description of such agglomeration processes should use local struc- 
tural configurations (LSC) defined by the concentration, corresponding to short-range order 
(SRO), and consequently to a random network description of the glass. We should stress 
here that the next level of description, using intermediate range order, is very similar to 
the SRO constructioni. The LSC's in GcxSci-x binary can for instance be the germanium 
tetrahedrally coordinated to selenium atoms, or in silica based glasses the silicon tetrahedra 
Q'-'^^ (the subscript k refers to the number of bonding oxygens on each S'z04/2 tetrahedron). 

Consider a typical cluster with a certain LSC distribution {p^i}i=i..N- As long as the 
viscosity is not too high, other LSC can stick on this typical cluster, creating new covalent 
bonds i-j with probability: 



where Wij is a statistical factor corresponding to the number of equivalent ways to stick a 
LSC z on a LSC j being part of a cluster {Wij is thus related to the coordination numbers 
TTii and TTij of the LSC) and Eij is the i-j LSC bond energy . Z normalizes the bond 
distribution. The creation of these new bonds produces a local variation in the probability 
(or concentration) distribution of the cluster and it can be encoded in the following master 
equation: 
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where r represents the mean agglomeration time and ( |2.2| ) represents a system of (A^ — 1) 
non hnear differential equations. At solidification temperature Tg {Tg = for a crystal or 
Tg = Tg for a glass) one should reach a stationary state and dpi/dt = 0, i.e. the variation of 
local probability distribution should be minimized. 

How can we distinguish a glass from a crystal ? Imagine that a local fluctuation appears 
in the vicinity of a stationary solution, satisfying a linearized version of equ. (|2.2|) . As usual, 
we can distinguish three types of singular points by means of the linearization. If all the 
roots of the characteristic equation of the linearized system have a negative real part, the 
solution is a stable attractor, i.e. it will show the preferential agglomeration process, which 
corresponds to nucleation of the crystal. Thus there will be no possibility for a fluctuation 
to grow, and we identify the stable stationary solution with a crystal, and T, = T^. If all the 
real parts are positive, one gets an unstable stationary solution. If both are present, a saddle 
point solution is obtained (flg. 1). The glass correponds to the latter characteristic, because 
it is neither a stable nor an unstable system, it has metastable character, and Tg = Tg. 
There is indeed still a chance for the system to escape (in other words for a fluctuation to 
grow) from the stationary saddle point and to fall on the stable crystalline attractor, which 
always happen experimentally when a glass is annealed. 

One can extend such a description to binary and ternary glasses as well. First, we 
consider the case when N = 2, i.e. when there are only two different types of LSC We 
denote them by A and B with their respective coordination numbers and m^. The 
system ( p.2|) reduces then to a single equation with only one variable, p° = x, e.g. the 
probability of occurence of the LSC B (and set equal to the concentration of B species). 
The solution of (|2.2| ) yields: 
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(2.3) 



where = Eab — Eaa is introduced when computing the probabilities paa and pab from 
equ. ( p^.l| ) (we have neglected the possibility of BB bonds because we are dealing in the fol- 
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lowing with low modified glasses only, and below the stoichiometric composition). One can 
see in the second part of equation ( |2.3| ) that the relationship can be made parameter- free, 
by considering the limit x ~ 0, when Tg ~ Tq. The initial network (with glass transition 
temperature Tq) is made of A LSC only (e.g. the selenium network in GcxSei-x systems) 
and one gets from the first part of equ. (PTB|): = fc^To In [m^/m^]. In order to predict 
the glass transition temperature in a binary glass Ai^^Bx, there is just need of the coor- 
dination number of the involved LSC's and the initial glass transition temperature Tq of 
the A network. Finally, we can obtain a linear equation at the very beginning of structural 
modification: 

(2.4) 
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The application to ternary glass networks Ai^^-yBxCy is slightly different, because when 
= 3, there are two non-linear equations to solve in terms of two probabilities p% = x 
and Pq = y. However, one obtains a saddle point solution from (|2.2D, yielding again a 
parameter-free relationship between x, y and Tg, because the new bond energy differences 
Ac = Eac-Eaa = kBTolnlmc/mA] and Ebc — Eaa = ^bTJ) In [mcms/m^] are determined 
again from boundary conditions (from the binary AC glass for the former, similarly to A^, 
from the binary slope equation for the latter)!. One interesting quantity in such systems 
(and in multicomponent chalcogenides) is the average coordination number, defined by f = 
771^(1 — X — y) + uibx + mcy (and rriA = 2). From the saddle point solution of equ. (p. 21), 
we obtain a relationship between f and Tg, to be compared with experiment: 



2mBTnc 



rriBmca'y^'ya — 7 — a) + 2rccP'{^ ~ 7) + '^'f'Bl'^i)- — o) 



(2.5) 

2rca — 2rB^ — tbTcoi'^ ) + SrBrcce'y 

where a = (2/mc)^'^°^'^^^ and 7 = (2/mB)^'^°^'^^^ ■ The slope in the limit f = 2 (i.e. x=y=0) 
has also a simple expression: 
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III. COMPARISON WITH EXPERIMENTAL DATA 



The obtained relationships ( p.4|) and ( |2.5| ) can be compared to the experimentally mea- 
sured glass transition temperatures in binary, ternary and multicomponent chalcogenide 
glass systems. 

Given the initial glass transition temperature To of vitreous sulphur (245 iS), selenium 
(316 K^) and tellurium (343 K extrapolated from the data ini), we have plotted the equa- 
tions ( ^.41) and ( ^.51 ) for chalcogenides including elements of Group IV and V. We can see 
that equation ( ^.41 ) predicts the Tg{x) trend at low concentration for all the binary systems 
IV- VI and V-VI systems displayed (fig. 2). From obvious structural considerations, we can 
insert in equ. (|2^ ) the value tjib = 4 (Group IV) or uib = 5 (Group V), = 2 and To, to 
be compared with the plotted experimental measurements on glass transition temperatures. 
The prediction gives also an indirect evidence of the stiffness transition (occuring at f = 2.4 
following the theory of Phillips and Thorpe!) . For f > 2.4 {x > 0.2 in IV- VI glasses), the 
network looses its random character and chemical ordering occurs, due to the chemical sta- 
bility composition at x=0.333. Thus the description in terms of a random network of A — A 
and A — B bonds should fail at this concentration. This is clearly seen for the GcxSci-x 



data, which start to deviate from the equation ( |2.4| ) at x = 0.18, and even more at x = 0.24, 
consistently with Mossbauer spectroscopy^. The sulfide system behaves very similarly, as 
seen on fig. 2 (and still tb = 4), although the structure of the initial glass (x=0) is rather 
different (chains and 5*8 rings). The addition of germanium leads to a random network com- 
posed of GeS'4/2 linkages between S chains and rings. Note that there is no deviation for the 
Ge-Te compound at x=0.2. This can be related to the fact that c — GeTe2 does not exist 
(in contrast with the existence of c — SiSe2, c — GeSe2, etc.) and probably that chemical 
ordering probably does not occur at the same concentration. As a consequence, the network 
of GcxTei-x can be thought as random. 

Agreement of the prediction with experimental measurements is also obtained in V-VI 



network glasses. For all the systems displayed in figure 2, the slope equation (|2.4|) gives the 
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correct trend in the variation of the glass transition temperature with network modification. 



Deviation of the stochastic prediction of ( |2.4| ) is here also supposed to occur at f = 2.4, 
corresponding to a; = 0.4 in the As-Se system. There is some evidence that the Bismuth 
atom could be five-fold coordinated in binary Bi-Se glassed. If this is the case, the stiffness 
transition may occur x ~ 0.13 (corresponding to f = 2.4 and may be observable in the Tg 
data, as well as the deviation of the stochastic prediction of equ. ( p.4D ). We can observe the 
same kind of agreement in ternary chalcogenides, where the glass transition temperature is 
given as a function of the average coordination number by equ. ( p.5|) or ( p.6|) . The figure 3 
shows that in very different systems (germanium chalcohalide glasses, thus = 4, mc = 1, 
and telluride glasses) the network is random for f < 2.4 and the description in terms of A-A, 
A-B, A-C and B-C bonds only is accurate. The same deviation is observable at f = 2.4 
which can again be interpreted by the occurence of chemical ordering due to the stiffness 
transition (e.g. occurence of a Sh2Te^ phase in the Ge-Sb-Te system). However, we should 
point out that experimental measurements in the low modified region of telluride systems 
(lower plot fig. 3) should be realized in order to definitely confirm the prediction. Other 
systems exhibit the same universal trends in the glass transition temperature variation and 
the construction can be extended to quaternary and multicomponent chalcogenide systems 
in a simple fashion!. 
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FIG. 1. The right-hand side of one of the equation ( p. 21) for = 3 in a simple polygon moc 
The system out of equilibrium can fall on the attractive = 1 stationary solution (crystallization, 
and Tg = Tm) and never on p2 = 1 or = 1. In some situations, the liquid can stay in the 
metastable state characterized by the saddle point solution at {pi = 0.6, p2 = 0.2, p^ = 0.2) and 
Ts = Tg. Note that the plot has to be truncated in order to have pi + P2 + = ^■ 
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FIG. 2. Binary IV- VI (upper plot) and V-VI ( lower plot) chalcogenide glasses (e.g. Ge^5ei- 



The lines represent the slope equation (2.4) with rriA = 2 and niB = 4 for the IV- VI glasses, and 
niB = 5 for V-VI glasses. Data have been displaced by 200 K and 100 K for a clearer presentationll. 
The vertical shaded line corresponds to the critical average coordination number fc = 2.4, predicted 
by Phillips and Thorpei. 
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FIG. 3. Glass transition temperature in ternary chalcogenides as a function of average coordination 



number f. Upper plot: chalcohalide glasses. Lower plot: ternary tellurides. The curves represent equ. (2.1) 
with rriA ~ 2 and ms and mc inserted following the Group of the Periodic Table. Data have been displaced 
for simplicity and are taken fromS 
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